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Stabilizing a Molecular Switch at Solid Surfaces: 
A Density-Functional Theory Study of Azobenzene at Cu(lll), Ag(lll), and Au(lll) 
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We present a density-functional theory trend study addressing the binding of the trans-cis confor- 
mational switch azobenzene (C6H5-N=N-C6H5) at three coinage metal surfaces. From the reported 
detailed energetic, geometric, and electronic structure data we conclude that the governing factor 
for the molecule-surface interaction is a competition between covalent bonding of the central azo 
(-N=N-) bridge on the one hand and the surface interaction of the two closed-shell phenyl (-CeHs) 
rings on the other. With respect to this factor the cis conformer exhibits a more favorable gas-phase 
geometric structure and is thus more stabilized at the studied surfaces. With the overall binding 
still rather weak the relative stability of the two isomers is thereby reduced at Ag(lll) and Au(lll). 
This is significantly different at Cu(lll), where the cis bonding is strong enough to even reverse 
the gas-phase energetic order at the level of the employed semi-local electronic exchange and cor- 
relation (xc) functional. While this actual reversal may well be affected by the deficiencies due to 
the approximate xc treatment, we critically discuss that the rationalization of the general effect of 
the surface on the meta-stable molecular states is quite robust. This should equally hold for the 
presented analysis of recent tip-manipulation and photo-excitation isomerization experiments from 
the view point of the derived bonding mechanism. 

PACS numbers: 68.43.Bc,71.15.Mb 
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I. INTRODUCTION 

In view of the rapidly advancing miniaturization in mi- 
croelectronics and sensing, molecules are envisioned as 
fundamental building blocks in a future "molecular nan- 
otechnology" . Since controlled switching between defined 
states is a crucial basis component for storage and logic, 
molecules offering this functionality (e.g. through exter- 
nally induced changes between conformational isomers) 
attain a central importance. Considering contacting and 
defined integration into a larger framework, it is more 
precisely the function of the molecule when stabilized at a 
solid surface that is of key interest. While a large variety 
of molecules can be controUably switched in gas-phase or 
solution, still little is known about their function in such 
an adsorbed state. Suppression of the switching capa- 
bility e.g. due to steric hindrance is a possible scenario, 
but completely new properties under the influence of the 
solid surface are equally conceivable. The atomic-scale 
understanding necessary for a technological exploitation 
of such effects builds on a detailed structural and elec- 
tronic characterization of the adsorbed molecular switch, 
as well as its response to external stimuli like fields, forces 
or external currentsii 

As a prototypical conformational switch of modest 
complexity azobenzene (C6H5-N=N-C6H5) has been on 
the research agenda for many yearsjS, In solution it ex- 
hibits a reversible photo-isomerization between an ener- 
getically more stable planar trans geometry and a three- 
dimensional cis configuration, in which the planes of the 
two phenyl-rings are tilted with respect to each other. A 
common theme in existing experimental studies address- 
ing the function in surface mounted geometries is that a 
too strong substrate interaction is suspected to cause ad- 



verse effects on the switching efficiency. This has either 
led to the use of ligands like alkanethiol chains to decou- 
ple the azobenzene moiety in vertical geomctriesii or to 
a focus on unreactive materials like close-packed surfaces 
of coinage metals and there in particular on Au(lll). 
On the latter surface bare azobenzene could be success- 
fully isomerized using tip-manipulation techniques^, but 
even in this alleged weak physisorption limit switching 
with light could not be achieved^. Reasoned as an ultra- 
fast quenching of the electronic excitation azobenzene- 
derivatives that are further lifted up from the surface ap- 
pear as a viable alternative and reversible photomechani- 
cal switching has indeed been reported for the functional- 
ized molecule tetra-tert-butyl azobenzene (TBA)^'^. The 
intricate role played even then by the metallic substrate 
is nevertheless exemplified by the fact that exactly for the 
same azobenzene-derivative, TEA, light-induced switch- 
ing could not be achieved at Ag(III)4 

As already stated, the characterization of the stable (or 
meta-stable, long-lived) molecular states at the surface 
is a necessary prerequisite to a detailed understanding of 
these experimental findings and of the actual excitation 
mechanism. While already less ambitious than a com- 
prehensive first-principles treatment of the entire switch- 
ing process, corresponding static electronic structure the- 
ory calculations still pose a considerable challenge. The 
large system sizes resulting from the sheer extension 
of the azobenzene molecule and the periodic supercell 
geometries dictated by the metallic bandstructure, can 
at present only be tackled by density-functional theory 
(DFT) with local- or semi-local exchange and correlation 
(xc) functionals. While this is the state of the art and 
the lowest level of theory at which one can still at least 
hope for the aspired quantitative and predictive-quality 
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modeling, it is clear from the start that the specificities of 
the molecule-surface binding directly challenge two well- 
known deficiencies of the named functionals: the spuri- 
ous self-interaction and the lack of long-range dispersive 
(van der Waals) attraction. With several frontier orbitals 
of different symmetry [vide infra) presumably involved 
in the bonding, the detrimental effect of self-interaction 
blurred orbital energies on fundamental quantities like 
the preferred adsorption site has already been demon- 
strated for much simpler adsorbates^iiSiiiii^ii^. Similarly 
well known is the important role played by long-range 
van der Waals interactions when aromatic 7r-like orbitals 
participate in the molecule-surface interactio n^^i^^ ?^^. 

The present trend study addressing the binding of 
azobenzene at the three coinage metal surfaces, Cu(lll), 
Ag(lll) and Au(lll) within large-scale DFT calculations 
correspondingly has a twofold focus. On the one hand, 
we critically discuss the obtained detailed energetic, geo- 
metric, and electronic structure data in the context of the 
sketched limitations of the employed semi-local xc func- 
tional, namely the generalized gradient approximation 
(GGA) functional due to Perdew, Burke and Ernzerhof 
(PBE)^^. With respect to the lack of long-range disper- 
sive interactions this includes the main conclusions from 
an assessment of this issue on the level of semi-empirical 
van der Waals correction scheme a^^'^^i^° , the details of 
which will be published in a consecutive paper—. On 
the other hand, we carve out some governing factors in 
the molecule-surface interaction that appear quite robust 
with respect to the employed approximate xc treatment. 
This concerns notably a rationalization of the obtained 
energetic ordering of the molecular states at the surface, 
and an analysis of the recent tip-manipulation and photo- 
excitation experiments from the view point of the derived 
bonding mechanism. As such, the reported data also pro- 
vides an important reference for future studies, in which 
we will investigate the differences induced by the engi- 
neered functionalization of this molecule to improve the 
switching efficiency. 



II. THEORY 

The DFT GGA-PBE calculations were carried out 
within the plane-wave implementation and using the 
standard library ultrasoft pseudopotentials-- as provided 
by the CASTEP code^^. Gas phase molecular geome- 
tries and properties were determined by optimizing iso- 
lated molecules in a rectangular (30 A x 40 A x 30 A) 
supercell, with F-point sampling and a plane-wave ki- 
netic energy cutoff of 800 eV. These supercell molecular 
calculations were complemented by computations of truly 
isolated molecules using the NWChem package^^ and the 
def-TZVP— basis set. While the geometries, orbital en- 
ergies and relative stabilities of neutral azobenzene ob- 
tained with both codes for GGA-PBE are identical on 
the accuracy level discussed here, we additionally used 
the molecular NWChem code for reference calculations 




FIG. 1: Schematic top view indicating the initial lateral posi- 
tion of the two azo-bridge N atoms in the six geometry opti- 
mizations carried out to determine the most stable adsorption 
site at each surface (see text for nomenclature and details) . N 
atoms, first- and second layer metal atoms are shown as black 
small circles, large light gray and large dark gray circles, re- 
spectively. 

using the hybrid B3LYP functional^^, as well as for the 
calculation of ionized molecular properties, cf. Table IIVI 
below. 

The molecule-surface system was modeled in supercell 
geometries, employing the GGA-PBE optimized lattice 
constants (ao(Cu) = 3.63 A, ao(Ag) = 4.14 A, ao(Au) = 
4.19 A) and, depending on the registry of the trans iso- 
mer with the substrate, large (6 x 3) or (6 x 4) (111) sur- 
face unit-cells to decouple the adsorbed molecule from 
its periodic images. No possibility for reconstruction 
was considered for both clean and adsorbate-covered sur- 
faces. The search for the most stable adsorption site at 
concomitant optimized adsorption geometries was per- 
formed using three layer slabs with adsorption on one 
side and allowing full relaxation of the topmost layer. 
Structure optimization relied on the BFGS algorithni^l 
as implemented in CASTEP, and the default convergence 
threshold on the maximum absolute ionic force compo- 
nent of 0.05 eV/A. The plane-wave kinetic energy cutoff 
was 350 eV, the vacuum region exceeded 13 A, and re- 
ciprocal space integration was done using (2x4x1) 
Monkhorst-Pack grids^®. Centrally targeted in these cal- 
culations was the relative energetic stability, which we 
found to be converged to ±20meV at these computa- 
tional settings. In order to determine the preferred ad- 
sorption site, both azobenzene isomers were optimized in 
six different adsorption geometries on all three surfaces 
as sketched in Fig. [1] The molecule was centered on 
the bridge (B), top (T) and hollow (H) surface sites, and 
shifted such that either the azo-bridge (— N — N— ) cen- 
ter of mass (designated as 1) or one of the azo-bridge N 
atoms (designated as 2) aligned with the adsorption site. 
In most optimizations, the molecule relaxed towards the 
1:1 metal-N atom coordinated Bl geometry (shown cen- 
tered at the bottom of Fig. [T]) , which at all three surfaces 
and for both isomers was also the energetically most fa- 
vorable among those found. 

For the thus defined most stable geometry the absolute 
adsorption energy, work function and projected density of 
states (PDOS) were determined at refined computational 
settings and a (6x3) (111) surface unit-cell. To also avoid 
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a finite dipole across the vacuum region we switched to 
thicker, inversion-symmetric seven layer slabs with ad- 
sorption on both sides. The vacuum was increased to 
22 A and the k-point density doubled. The adsorption 
energy is defined as 

£^ads = 2 [^azo@(lll) " '^'(111) ~ £^azo(gas)] , (1) 

where -Bazo@(iii) is the total energy of the relaxed 
azobenzene-surface system, £'(111) the total energy of the 
clean slab, -Eazo(gas) the total energy of the correspond- 
ing relaxed gas-phase isomer (all three computed at the 
same plane- wave cutoff), and the factor 1/2 accounts for 
the fact that adsorption is at both sides of the slab. The 
adsorption energy of either cis or trans isomer at the sur- 
face is thus measured relative to its stability in the gas- 
phase, and a negative sign indicates that adsorption is 
exothermic. Systematic convergence tests indicate that 
this central quantity is converged to within ±30 meV at 
the chosen settings. 

The work function of a particular surface structure is 
determined as the difference of the Fermi level and the 
value of the electrostatic potential at the center of the 
vacuum between two consecutive slabs, and we also ver- 
ified that this quantity is converged to within ±30meV 
at the chosen settings. As a technical point, we note that 
this required increasing the size of the grid on which the 
charge density is represented from the 1.8 kcut otherwise 
used, to the maximum of 2.0 kcut- Additionally, we had 
to lower the SCF convergence energy threshold from oth- 
erwise 10"'^eV/atom to 5 • lO^^eV/atom. 

The molecular PDOS^^ Pa{E) is defined by projec- 
tions of the Kohn-Sham (KS) states of the (adsorbed) 
molecule-surface system \(j)b,k) with KS Eigenvalues eh^k 
onto a certain KS orbital \ipa) of the isolated molecule in 
the adsorption geometry 

Pa{E) = Y,'^k\{Va\<t)Kk)? 5{E ~ebM) ■ (2) 

b,k 

The summation over band indices h and k-point indices 
k is weighted by the k-point-weights Wk so that an inte- 
gration of Pa[E) over E gives \{Lpa \ ^a)\^ if sufficiently 
many bands of the molecule-surface system are consid- 
ered. Conveniently, we used molecular KS orbitals |</5o,fc) 
which were calculated in the same supercell and with the 
same energy cut-off and k-point grid as the molecule- 
surface system to evaluate Eq. We verified that the 
artificial dispersion of these molecular KS orbitals is in- 
deed negligible within the desired accuracy of our compu- 
tational setup. In addition, Eq. Q was evaluated using 
the scalar product induced by the overlap matrix S due to 
the ultrasoft pseudopotentials^^. In the plots below, the 
(5-functions have finally been convoluted with a Gaussian 
function of 0.2 eV full width at half maximum for better 
visualization. 



TABLE I: Selected structural parameters of gas-phase 
azobenzene (AB) as defined in the text and in Fig. O be- 
low. Compared are the optimized values from the GGA-PBE 
calculations against higher-level theory (B3LYP (this work) 
and MP2/cc-pVTZ^) and experimen t ' . 
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III. RESULTS AND DISCUSSION 
A. Gas-phase azobenzene 

We begin our analysis with a short compilation of data 
for gas-phase azobenzene that is of relevance for the en- 
suing study of the surface mounted switch. With the 
molecular structure well summarized as two phenyl rings 
connected by an azo (— N = N— ) bridge, central geometric 
parameters are the azo-bridge bond length dNNj as well 
as the two angles that describe the orientation of each 
symmetry-equivalent phenyl-ring to the azo-bridge and 
of the phenyl-rings to each other, inversion angle a and 
dihedral angle w, respectively. Figure [3] below illustrates 
these parameters in the context of the surface adsorbed 
geometry, while Table H] compares their optimized gas- 
phase values obtained from the GGA-PBE calculations 
with higher-level theory and experiment. With the pos- 
sible exception of a slightly elongated azo-bridge bond 
length it is fair to say that with respect to the geomet- 
ric structure the employed GGA functional yields overall 
results that are as good as that of the hybrid B3LYP 
functional or the Hartree-Fock based M0ller-Plesset per- 
turbation theory (MP2). This favorable performance also 
carries over to the relative energetic stability of the two 
isomers. Here, we find the trans form to be more stable 
by 0.57 eV than the cis form, in good agreement to the 
values of 0.68 eV and ^ 0.6 eV from the B3LYP calcula- 
tions and from experiment^, respectively. 

As anticipated the description of the molecular fron- 
tier orbitals is much more problematic. Figure [2] shows 
representations of those KS orbitals that are of particu- 
lar interest both with respect to the bonding mechanism 
to the surface, and with respect to possible isomerization 
pathways. These range from the second highest occupied 
molecular orbital (HOMO-l), over HOMO and lowest un- 
occupied molecular orbital (LUMO) to LUMO-l-1. Com- 
pared to e.g. the much less self-interaction hampered 
B3LYP functional the GGA-PBE functional stiU yields 
the correct molecular character of these orbitals. In the 
higher-symmetry trans state the HOMO-l and LUMO 
thus contain large contributions of the tt and tt* orbital 
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FIG. 2: (Color online) GGA-PBE Kohn-Sham eigenvalues 
with respect to the vacuum level of the four frontier orbitals 
HOMO-1 (tt), homo (n), LUMO (tt*) and LUMO+1 (tt*) of 
gas-phase trans (solid lines) and cis (dashed lines) azobenzene. 
Additionally shown are corresponding orbital representations 
superimposed on the skeleton geometric structure (left: trans, 
right: cis). 



of the — N = N— moiety, respectively. The HOMO con- 
sists predominantly of the N lone pairs and the LUMO+1 
is centered on the tt* system of the phenyl rings. This 
character of the four frontier orbitals is also largely re- 
tained in the cis form, which is why they are often still 
classified as tt, n, or tt* despite the lack of symmetry of 
this isomer. While the ordering of these frontier orbitals 
is thus still identical in B3LYP and GGA-PBE calcula- 
tions, large differences exist in their actual energetic po- 
sition. This is dramatically illustrated by the computed 
HOMO-LUMO gap, which in GGA-PBE comes out at 
respectively 1.95 eV and 1.91 eV for trans and cis azoben- 
zene, and which is therewith about 2 eV smaller than at 
the B3LYP level (trans: 3.92 eV, cis: 3.76 eV). The cor- 
responding self-interaction driven underestimation of the 
HOMO-LUMO gap is typical for a GGA functionaP^ and 
the crucial question is how much of this error is retained 
in the adsorbed state. 



B. Adsorption geometry and energetics 

As already stated in Section II the most favorable ad- 
sorption site obtained at all three coinage metal surfaces 
corresponds to the — N=N— moiety bridging between two 
surface atoms as sketched in the Bl geometry of Fig. 
[TJ While in this adsorbed state the azo-bridge could in 
principle be vertically tilted with a concomitant symme- 
try breaking of the molecular structure as expressed by 
two different inversion angles a of the two phenyl-rings, 
we find that this is not the case. At all three surfaces, 
the azo-bridge lies flat at a height z above the surface, 
and the two phenyl-rings show identical inversion angles 
that are furthermore essentially unchanged compared to 
the gas-phase value, i.e. 115° and 124° for trans and cis 
azobenzene respectively, cf. Table ID Figure [3] illustrates 
the resulting adsorption geometry and the main struc- 




FIG. 3: (Color online) Perspective view of adsorbed trans and 
cis azobenzene in the most stable Bl adsorption site, cf. Fig. 
[T] illustrating the main parameters defining the adsorption 
geometry (see text): The azo-bridge bond length dNN and 
the vertical height z of the azo-bridge above the surface plane 
defined by the average height of the first-layer atoms. Fur- 
thermore the azo-bridge to phenyl-ring bond angles a and the 
azo-bridge dihedral angle uj. The latter is defined as the small- 
est angle formed between the — C— N— bond to one phenyl- 
ring and the plane spanned by the — N = N— bridge and the 
— N— C— bond to the other phenyl-ring. 



TABLE II: Structural parameters as defined in Fig. O for 
trans and cis azobenzene (AB) in the relaxed most stable 
adsorption geometry at the three coinage metal surfaces and 
compared to the optimized gas-phase structure. 
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3.50 


1.30 
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1.29 


18 



tural parameters discussed already in the context of the 
gas-phase molecule. 

If we start the detailed discussion with trans azoben- 
zene at Au(lll) and Ag(lll), one can summarize the 
calculated geometry values listed in Table |TT] in that the 
molecule remains essentially flat and floats at a rather 
large height > 3.5 A above the surface. Essentially flat in 
this context means that the dihedral angle to remains in 
both cases zero as in the gas-phase. Due to some torsional 
twisting of the phenyl-rings, the variation of z-heights ex- 
hibited by the individual atoms within the molecule nev- 
ertheless amounts to Az — 0.34 A and 0.37 A on Au(lll) 
and Ag(lll), respectively. While this is relatively small, 
it could nevertheless be an important feature consider- 
ing that recent conformational dynamics calculations^ 
for a free azobenzene derivative found that a non-planar 
geometry is a necessary prerequisite to efficiently couple 
electronic excitations with intramolecular rotations. 

The comparatively weak bonding suggested by the 
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TABLE III; Adsorption energies iJads as defined in Eq. ([1} for 
trans and cis azobenzene at tfie three coinage metai surfaces, 
as well as the relative energetic stability of the two adsorbed 
isomers Ai? with Ai? > indicating a higher stability of the 
trans isomer. All values in eV. 





Sads (trans) 


-Eads(cis) 


AE 


Cu(lll) 


-0.27 


-1.08 


-0.24 


Ag(lll) 


-0.11 


-0.42 


0.26 


Au(lll) 


-0.12 


-0.27 


0.42 


gas-phase 






0.57 



large 2:-height is indeed confirmed by the actually calcu- 
lated exothermic, but small adsorption energy of about 
-0.1 eV at both surfaces, cf. Table Hill and is furthermore 
consistent with the virtually unchanged azo-bridge bond 
length of 1.30 A in the adsorption geometry, cf. Table HTl 
The correspondingly deduced weak physisorptive bond- 
ing scenario is, however, particularly prone to the promi- 
nent lack of dispersive interactions in the employed GGA- 
PBE functional as discussed in the introduction. With- 
out a corresponding attractive contribution, the z-height 
of the molecule is exclusively governed by the steep Pauli 
repulsion wall encountered by the phenyl-rings, prevent- 
ing a closer approach and therewith stronger bonding 
of the molecule to the surface. Without further quali- 
fication of such interactions we would therefore refrain 
from putting too much emphasis on the actually calcu- 
lated value of the z-height and the adsorption energy. 
In fact, our investigation addressing this issue at least 
on the level of semi-empirical van der Waals correction 
schemes indicates that the dispersive attraction induced 
lowering of the z-height could amount to up to ^ 0.7 A, 
with a concomitant significant increase in the adsorption 
energji2i. 

Nevertheless, even at such a reduced height, the planar 
trans geometry fixes the azo-bridge atoms at a distance 
to the surface that is outside the range that would enable 
the formation of stronger covalent bonds to the top-layer 
metal atoms, and is thereby the main limiting factor in 
the molecule-surface interaction. This problem does not 
exist for the cis isomer, where due to the bent molecular 
geometry the azo-bridge can come much closer to the sur- 
face with the phenyl-rings pointing away from the surface 
and thus at heights that are not yet in conflict with the 
repulsive wall, cf. Fig. [S] This picture is nicely reflected 
by the computed much smaller z-heights slightly in excess 
of 2 A, cf. Table HB that bring the azo-bridge atoms into 
bond distances to the directly coordinated surface metal 
atoms that are highly reminiscent of corresponding dis- 
tances in transition metal complexes^. As a consequence 
the actual intramolecular — N=N— bond gets activated as 
reflected by the computed slightly elongated bond length 
compiled in Table |TT1 In detail, the resulting optimum 
z-heights still seem to require some additional bending 
of the molecule to lift the phenyl-rings even further away 
from the surface as expressed by the somewhat increased 



dihedral angle uj for cis azobenzene at both surfaces, cf. 
Table HT] The cost for this bending partly compensates 
the gain from the formed azo-bridge surface bonds, but 
as apparent from Table IIIII the net resulting adsorption 
energy for the cis isomer is still larger than in the trans 
case, and correspondingly, we arrive at a decrease of the 
relative energetic stability of the two isomers. 

Nevertheless, at both Au(lll) and Ag(lll) the trans 
form is still the more stable one, consistent with the 
conclusion from several STM studies which interpreted 
the measured "dumbbell" appearance of the thermally 
more stable adsorbed molecular conflguration in terms 
of the trans isomer— i^i^Z. While DFT calculations based 
on a local-density approximation (LDA)'^- xc functional 
supporting one of these STM investigations^ also found 
the trans isomer to be more stable, it is important to 
note that their reported bond strength to the surface is 
qualitatively different to our findings. For both isomers 
at Au(lll) strong chemisorptive binding of more than 
1.5 eV was computed, with the most favorable cis adsorp- 
tion geometry in fact featuring one phenyl-ring parallel to 
the surface and one phenyl-ring standing upright.^ Sus- 
pecting this difference to predominantly arise out of the 
use of the different, yet equally standard xc functional, 
this underscores the initially made statements concerning 
the sensitivity of the description of the molecule-surface 
interaction to the approximate treatment of electronic ex- 
change and correlation in present-day functionals avail- 
able for huge supercell calculations. When accounting 
for the lacking dispersive attraction in the here employed 
GGA-PBE functional with semi-empirical van der Waals 
schemes, we obtain a maximum correction of the adsorp- 
tion energy of the order of 1 eV^^. While this brings the 
absolute value closer to the reported LDA result, the one 
phenyl-ring upright cis geometry is never found to be sta- 
ble, let alone that it is clear that the LDA gives the bond 
strength for the wrong reasons. In this respect, we are 
confident that the approach of dispersion corrected GGA- 
PBE energetics employed in this study provides a much 
more adequate description of the azobenzene bonding to 
coinage metal surfaces, and we will critically discuss this 
point in more detail in our consecutive publication^!. 

Moving over to Cu(lll) as the remaining of the three 
surfaces we find a much stronger interaction even at the 
GGA-PBE level. In contrast to the situation at Ag(lll) 
and Au(lll), both the trans and cis form of azobenzene 
adsorbs at a vertical z-height that brings the azo-bridge 
atoms to typical covalent N-Cu bond distances known 
from transition metal complexes^. While due to the 
smaller Cu radius the Pauli repulsion for the phenyl- 
rings sets in at a closer distance^^, such low z-heights are 
nevertheless incompatible with a fully planar trans ad- 
sorption geometry. As shown in Table |TT] this leads to a 
pronounced buckling of the molecule with a final dihedral 
angle lo that is almost equal to the one of the adsorbed 
cis isomer, with even the latter angle being significantly 
increased with respect to its gas-phase value. The cost of 
this severe distortion of trans azobenzene again directly 
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carries over to the final adsorption energy. While the 
strong surface-interaction indicated by the much elon- 
gated -N=N- moiety is directly visible in the strong i^ads 
of the cis isomer, the adsorption energy in the trans case 
is correspondingly much smaller despite an even more 
activated azo bond length. 

At Cu(lll) this difference in the net adsorption energy 
finally leads to a reversal of the relative stability, i.e. we 
obtain the cis isomer to be more stable by — 0.24eV than 
the trans configuration, which as such is an intriguing 
finding. It is again intuitively rationalized by the differ- 
ent preferred adsorption heights of the central azo-bridge 
and the closed-shell phenyl-rings, which can be much bet- 
ter acommodated in the three-dimensional cis form. For 
any strong enough azo-bridge to surface bonding this will 
naturally offset the higher gas-phase stability of the trans 
isomer, and as indicated by our calculations this seems to 
be the case at Cu(lll). Nevertheless, the obtained pref- 
erence by — 0.24eV is a small number and has to be seen 
in the perspective that the GGA-PBE functional does 
not account for possibly different van der Waals contri- 
butions to the adsorption energy of the two isomers. This 
is particularly critical, as it is clear that these contribu- 
tions will yield a larger stabilization for the more planar 
trans configuration, and would therewith counteract the 
covalent bond driven stability reversal. Our analysis of 
this issue on the level of the semi-empirical van der Waals 
correction schemes in fact indicates that the stability re- 
versal does not prevaii^i. Nevertheless, in view of the 
unsure applicability of these schemes to the adsorption 
at metal surfaces it remains for future higher-level the- 
ory studies or experiment to fully settle this point. Not 
withstanding, despite the also partly significant changes 
to the geometric structure induced by the semi-empirical 
van der Waals correction schemes, we verified that this 
has only a minor effect on the electronic structure dis- 
cussed in the next section. We correspondingly present 
this analysis at the optimized GGA-PBE geometries and 
only briefly comment on the changes obtained at the van 
der Waals corrected geometries. 



C. Electronic structure 

Important information concerning the bonding mecha- 
nism at the surface is already provided by the electroneg- 
ativity (EN) of the isolated molecule. For trans and 
cis azobenzene we compute the Mulliken EN, defined as 
the mean of electron affinity and ionization potential, as 
4.44 eV and 4.11 eV, respectively, cf. TableHVl If we start 
our discussion again with the adsorption at Au(lll) and 
correspondingly compare these values to the work func- 
tion of the clean surface, $au(iii) = 5-20 eV, this suggests 
some charge transfer from the molecule to the metal sub- 
strate, and more specifically slightly more charge transfer 
in case of the less electronegative cis isomer. The actu- 
ally calculated lowering of the work function for trans and 
cis azobenzene adsorption of — 0.37eV and — l.SOeV, re- 



TABLE IV: Electron affinity (EA), ionization potential (IP), 
and Mulliken electronegativity (EN) of the free trans and cis 
azobenzene molecule, compared to the work function of the 
clean surface ($), the change of the work function induced 
by the adsorbed molecule (A$), and the change of the work 
function after subtracting the change induced by the poten- 
tial drop through the molecular overlayer (A<1?*), see text. 
The last two values correspond to the molecular coverage as 
resulting from the employed surface unit cell. Nevertheless, 
since we use the same (6 x 3) unit cell for both isomers, the 
induced work function changes are directly comparable. All 
values are given in eV. 





trans azo 


cis azo 


EA 




1.06 




0.75 


IP 




7.82 




7.47 


EN 




4.44 




4.11 


Cu(lll) 


$ 




4.83 






A$ 


-0.69 




-0.96 




A$* 


-0.15 




+0.11 


Ag(lll) 


$ 




4.46 






A$ 


-0.26 




-0.77 




A$* 


-0.25 




+0.03 


Au(lll) 


$ 




5.20 






A$ 


-0.37 




-1.30 




A$* 


-0.36 




-0.55 



spectively, is consistent with this expectation, cf. Table 
IIVI Quantitatively, it is, however, important to realize 
that even without any charge transfer (or more generally 
charge rearrangement) upon adsorption a work function 
change can already result if the molecular overlayer as 
such exhibits a non-zero dipole moment^. In this re- 
spect, the two azobenzene isomers differ significantly in 
that the gas-phase trans configuration has no dipole mo- 
ment, whereas the cis configuration exhibits a finite cal- 
culated dipole moment of 3.23Debye. In order to dis- 
entangle the contributions to the work function change, 
we therefore compute the potential change across a free- 
standing molecular overlayer in exactly the same geo- 
metric structure as it has in the adsorbed state. Sub- 
tracting this term we arrive at the work function change 
A$* solely due to the charge rearrangement caused by 
the interaction of the overlayer with the metal substrate. 
As shown in Table llVi the values for A$^^^^^^^j are 
— 0.36eV and — 0.55eV for trans and cis azobenzene, re- 
spectively, which now fit almost too perfectly into the 
electronegativity picture of an only slightly larger charge 
transfer to the metal in case of the cis isomer. At this 
stage, we stress, however, the well-known limitations of 
the charge transfer notion at metallic surfaces.^^ Distin- 
guishing between an actual transfer of charge between 
molecule and substrate or a polarization of the molec- 
ular charge accompanied by a varying degree of charge 
displacement due to Pauli repulsion in the overlapping 
electron density tails of the bonding partners is to some 
degree semantics 42. We make no attempt to contribute 
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to this discussion here, and henceforth simply refer to 
charge transfer in quotes to emphasize its conceptual na- 
ture. 

Figure m analyses the azobenzene/Au(lll) interaction 
in terms of the density of states projected onto the molec- 
ular frontier orbitals. Similar to the discussion of the 
work function changes we aim to single out the effect 
of the charge rearrangement upon adsorption by also 
showing the level positions as they would result within a 
Schottky-Mott model, i.e. by simply aligning the vac- 
uum level of the free molecule and the clean surface. 
While this model in general ignores the formation of a 
surface dipole, we modify it for our purposes by including 
the above described potential change caused already by 
the dipole moment of the free-standing molecular over- 
layer. Any shift between the orbital energy in this modi- 
fied Schottky-Mott model and in the fully self-consistent 
calculation can this way be attributed to the charge re- 
arrangement caused by the molecule-surface interaction. 
Starting with the PDOS of trans azobenzene at Au(lll) 
shown in the left upper panel of Fig. 31 we see the four 
projected frontier orbitals, HOMO-1, HOMO, LUMO, 
and LUMO-f 1 broadened with hardly any degree of inter- 
mixing with Au d-states. Compared to their position in 
the modified Schottky-Mott model all levels are roughly 
equally shifted by about 0.8 eV to lower energies, which 
is the direction one would also expect from the afore dis- 
cussed "charge transfer" into the metal substrate, but 
which could equally well merely result from the lower- 
ing in the surface potential. Despite the down-shift, the 
LUMO still remains well above the Fermi-level. Both this 
and the actually resulting HOMO and LUMO positions 
at — 0.7eV and -1-1.2 eV respectively are in reasonable 
agreement with scanning-tunneling spectroscopy (STS) 
measurements by Choi et al^. This suggests that the 
sizable self-interaction induced error in the description 
of the HOMO-LUMO gap of the free molecule, cf. Sec- 
tion IIIA, is to some extent lifted, and that the orbital 
energies of the surface adsorbed molecule may be rather 
adequately described at the GGA-PBE level. 

In case of cis azobenzene adsorption at Au(lll) the 
actual level positions are in general at lower energies, cf. 
Fig. m This is primarily due to the fact that because of 
the finite potential change caused by the free-standing 
molecular overlayer already the level positions in the 
modified Schottky-Mott model lie at lower energies. 
With a comparable down-shift of the self-consistently cal- 
culated frontier orbitals as for trans azobenzene {vide 
infra), this has as one consequence that the projected 
LUMO lies now just above the Fermi-level. As a sec- 
ond consequence, both the HOMO and HOMO-1 level 
fall right into the energy range of the Au d-band with 
a concomitant stronger hybridization. The down-shift of 
the LUMO is with about 1.0 eV even slightly stronger 
than for the trans isomer, again in line with the above 
discussed picture of a slightly stronger "charge transfer" 
as inferred from electronegativity differences and work 
function changes, or alternatively with a stronger lower- 



ing in the surface potential due to the lower adsorption 
height of this isomer. The exception of the much smaller 
shift of the LUMO-l-1 level of only 0.4 eV is in this respect 
primarily due to the distorted adsorption geometry: The 
computed LUMO+1 for a free molecule in this distorted 
geometry is higher in energy than in the relaxed gas- 
phase geometry, and if one takes this into account, one 
arrives at an effective down-shift of the LUMO-l-1 that is 
comparable to that of the other frontier orbitals. 

Moving to the azobenzene adsorption at Ag(lll) we 
find essentially the same picture as obtained at Au(lll), 
yet with subtle and intriguing modifications. With a 
much lower work function of fE'Ag(iii) = 4.46 eV the 
electronegativity difference to the azobenzene molecule 
is much smaller, and one would correspondingly expect 
a smaller "charge transfer". Consistently, the reduced 
work function change A<I>* upon trans azobenzene ad- 
sorption, i.e. after correcting for the dipole moment of 
the free-standing molecular overlayer, is with — 0.23eV 
smaller than the corresponding value at Au(lll), cf. Ta- 
ble IIVI Also the again roughly equal down-shift of the 
projected frontier orbitals compared to their position in 
the modified Schottky-Mott model shown in Fig. 0] is 
with about 0.6 eV smaller than the corresponding value 
of 0.8 eV for trans adsorption at Au(lll). Neverthe- 
less, due to the much smaller Ag(lll) work function the 
LUMO still ends up closer to the Fermi-level than was 
the case at Au(lll), and is further lowered to just above 
E-p in the van der Waals corrected geometry. Neverthe- 
less, with the Ag d-band starting at lower energies not 
even the HOMO-1 has significant overlap with the metal 
d-states, and we find all four frontier orbitals largely re- 
taining their molecular character. While all of these find- 
ings fit perfectly into the overall picture developed at 
Au(lll), the deviations start when turning to the cis 
azobenzene adsorption. From the smaller electronega- 
tivity of this isomer, one would expect a slightly larger 
"charge transfer" than in the case of trans azobenzene, 
and correspondingly a somewhat larger reduced work 
function change A<i>* as was the case at Au(lll). How- 
ever, actually calculated after correcting for the change 
induced by the free-standing molecular overlayer is an 
essentially zero A$*, cf. Table HVl We attribute this to 
the closeness of the LUMO to the Fermi-level already in 
the modified Schottky-Mott model as apparent in Fig. |4l 
A "charge transfer" or surface potential induced down- 
shift of this orbital leads therefore to a significant popula- 
tion of this formerly unoccupied level with a concomitant 
charge back-donation into the molecule counteracting the 
work function lowering. 

This more complex interaction is finally also obtained 
for trans and cis azobenzene at Cu(lll), where we now 
find a significant LUMO population for both isomers. 
That this results also for the trans case is partly already 
due to the strongly distorted adsorption geometries de- 
scribed in Section IIIB above. These distorted geometries 
exhibit a much increased dipole moment, which offsets 
the somewhat larger Cu(lll) work function, cf. Table 
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FIG. 4: (Color online) Density of states (DOS) projected onto the frontier orbitals of gas-phase azobenzene: HOMO-1 (blue 
dotted line), HOMO (green dashed line), LUMO (red solid line), LUMO+l (magenta dash-dotted line), see Fig. [2] for a 
representation of these orbitals in the free molecule. Additionally shown is the total DOS of the employed supercell downscaled 
by a factor 75 (thin black line), as well as the position of the corresponding frontier orbitals (vertical sticks at the upper 
a;-axis) in a modified Schottky-Mott picture that neglects the formation of a surface dipole through charge rearrangement upon 
adsorption (see text). Upper panel: Au(lll). Middle panel: Ag(lll), Lower panel: Cu(lll). For all three surfaces adsorption 
of the trans isomer is shown on the left, and adsorption of the cis isomer on the right. 



mil and leads to level positions in the modified Schottky- 
Mott model that are as low as in the case of Ag(lll). 
With the higher-lying Cu d-band and the much shorter 
distance of the adsorbed molecule above the Cu(lll) sur- 
face this enables a much stronger hybridization of the 
frontier orbitals than at the other two coinage metal 
surfaces, which is even further increased at the van der 



Waals corrected geometry. As a result HOMO-1, HOMO 
and LUMO all form broad bands with mixed metallic 
and molecular character and the counteracting donation 
and back-donation contributions result in overall small 
reduced work function changes A$* for both isomers, cf. 
Table IVl 
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IV. CONCLUSIONS 

We have presented a DFT GGA-PBE trend study of 
the adsorption of azobenzene at the three coinage metal 
surfaces, Cu(lll), Ag(lll), and Au(lll), providing de- 
tailed energetic, geometric, and electronic structure data 
of the meta-stable surface molecular states. As a cen- 
tral insight the calculations identify the conflicting bond- 
ing interests of the azo-bridge moiety on the one side 
and of the phenyl-rings on the other side as the govern- 
ing factor for the molecule-surface interaction. While 
the phenyl-rings prefer a rather large, Pauli-repulsion 
dictated z-height above the surface, much shorter dis- 
tances are required to enable the formation of stronger 
covalent N-metal bonds. For the planar trans config- 
uration these two demands are incompatible and as a 
result we find the molecule either floating at a larger 
z-height above Au(lll) and Ag(lll), or severely buck- 
led and distorted at a shorter z-height above Cu(lll). 
The three-dimensional cis configuration does not suffer 
from this conflict and the actually calculated adsorption 
height corresponds at all three surfaces to N bond lengths 
to the directly coordinated surface metal atoms that are 
highly reminiscent of corresponding distances in transi- 
tion metal complexes. 

Due to the more favorable gas-phase molecular geom- 
etry adsorption leads thus naturally to a larger stabiliza- 
tion of the cis isomer. At the employed GGA-PBE level 
we correspondingly compute a diminished relative stabil- 
ity with the trans isomer still more stable at Au(lll) and 
Ag(lll), and an actual reversal of the gas-phase energetic 
order at Cu(lll). This is traced back to the more reac- 
tive Cu d-band in conjunction with the smaller Cu atomic 
radius, which together enable the formation of strong co- 
valent N-Cu bonds at a low adsorption height and a con- 
comitantly reduced Pauli-repulsion of the phenyl-rings. 
From a critical discussion of these findings including a 
semi-empirical account of long-range van der Waals forces 
we conclude that the stronger dispersive attraction in 
case of the planar trans isomer largely counteracts the 
covalent bond driven preferential stabilization of the cis 
isomer. Whether the stability reversal of cis and trans 
azobenzene obtained with GGA-PBE at Cu(lll) indeed 
prevails remains thus to be confirmed by higher-level the- 



ory or experiment. 

Comparison to scarce experimental data suggests that 
the energetic positions of the frontier orbitals in the sur- 
face mounted geometry are reasonably well described 
to also attempt some conclusions with respect to the 
bonding and isomerization mechanisms. Light-induced 
switching in gas-phase and solution is conventionally dis- 
cussed in terms of photo-excitation of the HOMO-1 to 
LUMO (tt ^ TT*) and HOMO to LUMO (n tt*) reso- 
nances in the trans and cis states, respectively.^ In this 
respect, the obtained noticeable broadening of these lev- 
els at all three surfaces casts severe doubts on the effi- 
ciency of these conventional mechanisms, and could ra- 
tionalize the unsuccessful attempts to photo-isomerize 
bare azobenzene at Au(lll)^. On the other hand, 
the actual degree of hybridization of the occupied or- 
bitals could just be the key to the function, considering 
the HOMO-involving substrate-mediated charge-transfer 
process suggested recently in the context of the reversible 
photo-induced isomerization of the azobenzene-derivative 
TBA at Au(lll}^. In this regard, we expect further in- 
sight from a detailed comparison of the present reference 
data for bare azobenzene with corresponding data for dif- 
ferently functionalized azobenzene molecules. A last in- 
triguing aspect is the partial occupation of the LUMO in 
both isomers at Cu(lll) and in the cis isomer at Ag(lll) 
with a corresponding back-donation contribution to the 
bonding. It is thus only at Au(lll) that the LUMO po- 
sition remains distinctly above the metal Fermi-level for 
both isomers. This has quite some bearing on the neg- 
ative ion resonance mechanism discussed in the context 
of tip-manipulation techniques^ and could explain why 
the hitherto only such isomerization has been reported 
for this surface^. 
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